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The scaled harmonic oscillator basis (SHO) is derived by a local scaling- 



' point transformation of the spherical harmonic oscillator radial wave func- 

tions. The unitary scaling transformation produces a basis with improved 
asymptotic properties. The SHO basis is employed in the solution of the rel- 
\ ativistic Hartree-Bogoliubov (RHB) equations in configurational space. The 

\ model is applied in the self-consistent mean-field approximation to the de- 

scription of the neutron halo in Ne isotopes. It is shown that an expansion 
of nucleon spinors and mean-field potentials in the SHO basis reproduces the 
pH \ asymptotic properties of neutron densities calculated by finite element dis- 

cretization in the coordinate space. In the RHB description of neutron skins 
5h ' and halos, SHO bases in two or three dimensions can be a useful alternative 

to technically complicated solutions on a mesh in coordinate space. 



I. INTRODUCTION 



Theoretical models based on the relativistic mean-field approximation have been success- 
fully applied in the description of a variety of nuclear structure phenomena in spherical and 
deformed /^-stable nuclei , and more recently in studies of exotic nuclei far from the valley 
of beta stability Short lived isotopes with extreme neutron - to - proton ratios display 

unique properties. On the neutron-rich side in particular, extremely diffuse and spatially 
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extended neutron densities are found in nuclei close to the particle drip lines. Densities with 
large spatial extensions result from the extremely weak binding of the outermost nucleons 
and the coupling between bound states and the particle continuum. The theoretical descrip- 
tion of the appearance of neutron skins, or eventually neutron halos, necessitates a careful 
treatment of the asymptotic part of the nucleonic densities. While calculations based on 
modern Monte Carlo shell-model techniques provide an accurate description of the struc- 
ture of relatively light exotic nuclei, for medium-heavy and heavy systems the only viable 
approach at present are mean-field models. 

In addition to the self-consistent mean-field potentials, pairing correlations have to be 
taken into account in order to describe ground-state properties of open shell nuclei. For 
strongly bound systems pairing can be included in the form of a simple BCS approximation 
in the valence shell. For drip line nuclei on the other hand, the separation energy of the 
last nucleons can become extremely small. The Fermi level is found close to the particle 
continuum, and the lowest particle-hole or particle-particle modes are embedded in the con- 
tinuum. The BCS pairing model presents only a poor approximation and predicts unbound 
systems. A correct description of the scattering of nucleonic pairs from bound states to 
the positive energy continuum, and therefore of the formation of spatially extended nucle- 
onic densities, necessitates a unified treatment of mean-field and pairing correlations. In 
the non relativistic framework the structure of spherical exotic neutron-rich, as well as pro- 
ton drip-line nuclei, has been described in the Hartree-Fock-Bogoliubov (HFB) theory in 
coordinate space [§0. Among models based on quantum hadrodynamics, the relativistic 
Hartree-Bogoliubov (RHB) model |]3|JT0[1 has been very successfully employed in calculations 
of ground-state properties of spherical exotic nuclei on both sides of the valley of /^-stability. 
For /5-stable nuclei the traditional Hartree-Fock models and the relativistic mean-field the- 
ory predict very similar ground-state densities. On the other hand, recent applications to 
the structure of drip-line nuclei have shown significant differences, especially in the isospin 
dependence of the spin-orbit potential. 

In Ref. 1^ it was emphasized that for a correct theoretical description of nuclei with 
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large neutron excess, the model must include: (i) a self-consistent solution for the mean- 
field in order to obtain the correct increase of the surface diffuseness, (ii) the correct isospin 
dependence of the spin-orbit term in the mean-field potential, (iii) a unified self-consistent 
description of pairing correlations, (iv) coordinate space solutions for the equations that 
describe the coupling between bound and continuum states. However, the coordinate space 
representation of HFB or RHB equations still presents serious difficulties in the calculation 
of ground-state properties of deformed nuclei. For finite-range interactions the technical and 
numerical problems on a two-dimensional mesh are so involved that reliable self-consistent 
calculations in coordinate space should not be expected soon, and therefore one has to 
look for alternative temporary solutions. In the configurational representation the HFB or 
RHB systems of partial differential equations have been solved by expanding the nucleon 
wave functions and mean-field potentials in a large oscillator basis. In many applications 
an expansion of the wave functions in an appropriate harmonic oscillator (HO) basis of 
spherical or axial symmetry provides a satisfactory level of accuracy. For nuclei at the drip 
lines the expansion in the localized oscillator basis presents only a poor approximation for 
the continuum states. Oscillator expansions produce densities which decrease to steeply in 
the asymptotic region at large distance from the center of the nucleus. The calculated rms 
radii cannot reproduce the experimental values, especially for halo nuclei. 

In a recent work |[Tl[] a new configurational representation basis has been suggested, based 



on a unitary transformation of the spherical HO basis. The new basis is derived from the 



standard oscillator basis by performing a local-scaling point coordinate transformation [[12 
The precise form of the transformation is determined by the desired asymptotic behavior of 
the densities. At the same time it preserves many useful properties of the HO wave functions. 
Using the new basis, characteristics of weakly bound orbitals for a square-well potential 
were analyzed, and the ground-state properties of some spherical nuclei were calculated in 
the framework of the energy density functional approach. It was shown that properties of 
spatially extended states can be described using a rather modest number of basis states. 
The new basis, it was suggested, might present an interesting alternative to the algorithms 



that are being developed for the coordinate space solutions of the HFB equations. In the 
present work we investigate a particular scaled HO basis which might be appropriate for the 
solution of the coupled system of Dirac-Hartree-Bogoliubov integro-differential eigenvalue 
equations, and Klein-Gordon equations for the meson fields. 

In Sec. 2 we present an outline of the relativistic Hartree-Bogoliubov theory and discuss 
the various approximations. The local-scaling transformed spherical HO basis is introduced 
in Sec. 3, and results of calculations in the configurational representation for a series of 
Ne isotopes are compared with solutions obtained by discretization of the system of RHB 
equations on the finite element mesh in coordinate space PJl^. A summary of our results 
is presented in Sec. 4. 



II. RELATIVISTIC HARTREE-BOGOLIUBOV MODEL 

The relativistic Hartree-Bogoliubov theory provides a unified description of mean-field 
and pairing correlations in nuclei. In this section we briefly outline the essential characteris- 
tics of the model that will be used in the calculation of ground-state properties. A detailed 
description of the model can be found in Refs. 0,0]. The ground state of a nucleus 1$ > is 
described as vacuum with respect to independent quasi-particle operators. The generalized 
single-nucleon Hamiltonian contains two average potentials: the self-consistent mean-field F 
which encloses all the long range particle-hole (ph) correlations, and a pairing field A which 
sums up the particle-particle (pp) correlations. The variation of the energy functional with 
respect to the hermitian density matrix p and the antisymmetric pairing tensor k, produces 
the single quasi-particle eigenvalue equations. The model was formulated in Ref . [|10| . In the 



Hartree approximation for the self-consistent mean field, the relativistic Hartree-Bogoliubov 
(RHB) equations read 

hD-m-X A \ /f/fc(r)\ /f/fc(r)\ 

-A* -hD + m + x) \Vk{r)) \Vkir) ) ' 

where /id is the single-nucleon Dirac Hamiltonian and m is the nucleon mass 
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hj, = -ta-V + /3(m + g,a(r)) + g^r3^°(r) + g.p" (r) + e-^i^A°(r). (2) 

The chemical potential A has to be determined by the particle number subsidiary condition 
in order that the expectation value of the particle number operator in the ground state 
equals the number of nucleons. The column vectors denote the quasi-particle spinors and 
Ej; are the quasi-particle energies. The Dirac Hamiltonian contains the mean-field potentials 
of the isoscalar scalar cr-meson, the isoscalar vector cj-meson, the isovector vector p-meson, 
as well as the electrostatic potential. The RHB equations have to be solved self-consistently 
with potentials determined in the mean-field approximation from solutions of Klein-Gordon 
equations. The equation for the sigma meson contains the non-linear a self-interaction 



terms |Tj]. Because of charge conservation only the three-component of the isovector p- 
meson contributes. The source terms for the Klein-Gordon equations are calculated in the 
no -sea approximation. 

The pairing field A in ([^) reads 

A„b(r, r') = ^Y. ^abcd{r, r')Kcd(r, r'), (3) 

c,d 

where a, b, c, d denote quantum numbers that specify the single-nucleon states, Kbc(i(r, r') 
are matrix elements of a general two-body pairing interaction, and the pairing tensor is 
defined 

KUr,r')= Y: U:,{r)V,,ir'). (4) 

Ek>0 

The matrix elements in the pairing channel Vabcd can in principle be derived as a one- 
meson exchange interaction by eliminating the mesonic degrees of freedom in the model 
Lagrangian [|10|. However, the resulting pairing matrix elements are unrealistically large 
if for the coupling constants the standard parameter sets of the mean-field approximation 
are used. These parameters do not reproduce scattering data On the other hand, since we 
are using effective forces, there is no fundamental reason to have the same interaction both 
in the particle-hole and particle-particle channel. Although encouraging results have been 



reported in applications to nuclear matter, a microscopic and fully relativistic derivation 
of the pairing force starting from the Lagrangian of quantum hadrodynamics still cannot 
be applied to realistic nuclei. Therefore, it was suggested in Ref. [jl5[ that an appropriate 
approximation might be a phenomenological Gogny-type finite range interaction 



j=l,2 



-HiP^ -M.P^P^), (5) 



with the parameters /Xj, Wi, Bi, Hi and Mi {i = 1,2). The procedure requires no cut-off and 
provides a very reliable description of pairing properties in finite nuclei. 

The eigensolutions of Eq. (|l]) form a set of orthogonal and normalized single quasi- 
particle states. The corresponding eigenvalues are the single quasi-particle energies. The 
self-consistent iteration procedure is performed in the basis of quasi-particle states. The 
resulting quasi-particle eigenspectrum is then transformed into the canonical basis of single- 
particle states, in which the RHB ground-state takes the BCS form. The transformation 
determines the energies and occupation probabilities of the canonical states. 

The self-consistent solution of the Dirac-Hartree-Bogoliubov integro-differential eigen- 
value equations and Klein- Gordon equations for the meson fields determines the nuclear 
ground state. In Refs. [|2|-^ we used finite element methods in the coordinate space dis- 
cretization of the coupled system of equations. The model was applied to exotic neutron- 
rich nuclei. Very diffuse ground-state densities were described and the occurrence of neutron 
halo was investigated. In the description of not so exotic Ni and Sn isotopes 0, as well as 
for proton drip line nuclei , we followed the prescription of Ref. : the Dirac-Hartree- 



Bogoliubov equations and the equations for the meson fields were solved by expanding the 
nucleon spinors Uk{r) and Vfc(r), and the meson fields in a basis of spherical harmonic os- 
cillators for = 20 oscillator shells. For nuclei at the proton drip lines we verified the 
results by comparing calculations in configuration space with coordinate space solutions. 
Differences were negligible. As will be shown in the next section, things are very different 
on the neutron-rich side. 
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III. CALCULATIONS IN CONFIGURATIONAL REPRESENTATION 



A. The scaled harmonic oscillator basis 



In the construction of the local-scahng transformed harmonic oscillator basis we will use 
the local-scaling transformation method |T^. A local-scaling point coordinate transforma- 



tion (LST) is defined as 

r' = f(r)^f/(r). (6) 

The transformed radius vector has the same direction r = r/|r|, while its magnitude r' = /(r) 
depends on the scalar LST function, /(r) is assumed to be an increasing function of r, 
and /(0)=0. The set of invertible transformations of this type forms a LST group. The 
corresponding LST wave function can be expressed as 



^/(ri,r2, ...,rA) 



n 



i=i rf dn 



vI/(f(rO,f(r2),...,fM), (7) 



where \l/(ri, r2, ta) is a model A-particle wave function normalized to unity 

(^1^) = 1- (8) 

The transformation (|^) is unitary and therefore the LST wave function \&/(ri,r2, ...,r^) is 
also normalized to unity, independent of the choice of /(r). 

The local one-body density corresponding to a A-body wave function \1/ is 

p{r)=A J |^(r,r2,...,r^)pdr2,...,rfr^. (9) 

From the expression (0) it follows that there exists a simple relation between the local density 
P/(r) associated with the LST function and the density p(r) which corresponds to the 
prototypical model function 

Mr) - ^'-§^mr)). (10) 

This relation is particularly useful when the form of the density p/(r) is known, for example 
from experimental data. In that case Eq. ([TOf) becomes a first order nonlinear differential 



equation for the LST function /. For a system with spherical symmetry, p/, p, and / depend 
only on r=|r|, and Eq. (^) can be reduced to a nonlinear algebraic equation 

fir) 

J pf{u)u'^du = j p{u)v?du. (11) 



The solution is found subject to the boundary condition /(O) = 0. 

For shell-model or mean-field applications we have to consider the case when the model 
many-body wave function is a Slater determinant 

^(ri,r2,...,rA) = ^ det |0i(r,)|. (12) 

The single-particle functions 0i(r) form a complete set. The LST wave function is defined 
by the transformation (^, and of course it is written as a product state 

^/(ri,r2,...,r^) = ^ det |(/),(r,)|, (13) 



of the LST single-particle basis states 



(f(r)). (14) 



J.2 Qj. 

In the present work we consider systems with spherical symmetry. The single-particle 
basis reads 

(^i(r,(j,r) =R^^{r)Jji^{VL,a)xi/2mt(T)- (15) 
The angular part describes the coupling of orbital angular momentum and spin 

Ji{Vt,a) = Jjim^{Q,a) = {lmil/2ms\jmj)Yimi{^)Xi/2ms{<T)- (16) 

mims 

The harmonic oscillator radial functions 



are characterized by a single parameter: the oscillator length h = ^fi/muj. 

If applied in the description of very weakly bound nuclei, the main drawback of HO wave 
functions is that, because of their Gaussian asymptotics, they cannot describe the density 
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profiles far away from the center of the nucleus. In principle one could try to resolve this 
problem by including a very large number of oscillator shells, but extremely large bases are 
impossible to handle in more microscopic applications. Therefore, by means of eqs (0) and 
(p!5|) it might be useful to define a scaled harmonic oscillator (SHO) basis 



(f)i{r,a, t) =Rlf''{r)Jjim{n, a, r)xi/2m,(r). (18) 

The angular part of course does not change, while the radial functions R^^{r) are defined 
by the local-scaling transformation 



j,2 ^j. 



Knfir)). (19) 



The choice of the LST function 

r for r < a 

fir) = 



(20) 



1 /2 

a (8r /a — 8a/r + a^/r^ — 12 ln(r/a)) for r > a 

is motivated by the desired asymptotic form of the nucleonic densities. For this particular 
parametrization the derivatives of /(r) up to fourth order are continuous at the matching 
radius a, and /(r) is a monotonically increasing function of r. The transformation (|TU]) is 
unitary; the Jacobian ^-^^^^ is positive for all values of r. This makes the SHO basis com- 
plete, since it is unitary equivalent to the complete HO basis. The ansatz (^) for the LST 
function guarantees that all SHO basis states (|18D are localized in space. Asymptotically, 
the linear term in the expansion ( ]20|) dominates and /(r oo) ~ r^/^. Therefore, at large 
distances the radial functions decrease exponentially. The SHO functions are continuous 
(up to fourth derivatives) and depend on two positive parameters: the oscillator parameter 
h and the matching radius a. 

It should be emphasized that the simple parametrization (pO[) is motivated by our desire 



to keep the new basis as simple and practical as possible. In Ref. it has been shown 
that a three-parameter form of the LST function performs very well in actual calculations 
of ground-state properties of spherical nuclei. For light nuclei one of these parameters turns 



out to be very small. In the present application we choose it identically zero. The resulting 
matching conditions are satisfied analytically and the parametrization (^) is derived. Of 
course other parametrizations of /(r) are possible. There are many ways to optimize the 
LST function, depending on the actual physical problem under consideration. 

B. Neutron-rich Ne isotopes 

In Ref. [H we have applied the RHB model to the description of the neutron halo in the 
mass region above the s-d shell. Pairing correlations and the coupling to particle continuum 
states have been described by finite range two-body Gogny-type interaction. Finite element 
methods have been used in the coordinate space discretization of the coupled system of 
Dirac-Hartree-Bogoliubov and Klein-Gordon equations. As we have already emphasized, 
solutions in coordinate space appear to be essential for the correct description of the coupling 
between bound and continuum states. Using the parameter set NL3 [0 for the mean-field 
Lagrangian, and DIS parameters for the Gogny interaction, we have found evidence 
for the occurrence of multi-neutron halos in heavier Ne isotopes. We have shown that the 
properties of the lf-2p orbitals near the Fermi level and the neutron pairing interaction play 
a crucial role in the formation of the halo. In the present work we essentially repeat the 
calculations of Ref. , but in configurational representation with the SHO basis defined by 
the LS transformation ([l8|) . 

In Fig. ID we plot the proton and neutron rms radii for the Ne isotopes as functions of the 
mass number. Results of both calculations in the coordinate space and in the configurational 
representation (HO and SHO bases) are shown. Proton radii change very little and only 
display a slow increase. There is practically no difference between results of coordinate 

and configurational space calculations. Neutron radii of Ne isotopes follow the mean-field 

NV3 

curve up to N ^ 22. For larger values of N a sharp increase of neutron radii is 
observed (calculations in coordinate space). The matter radii of course follow this change 
and the spatial extension of heavier isotopes is substantially enlarged. The sudden increase 
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in neutron rms radii has been interpreted as evidence for the formation of a multi-particle 
halo 0. We notice that neutron radii calculated with the expansion in the HO basis cannot 
reproduce the coordinate space results. Results calculated with 12 and 20 oscillator shells are 
compared, and in both cases the radii follow the average mean-field curve. The neutron radii 
increase with the number of oscillator shells, but their values are still far from the results 
obtained in coordinate space. On the other hand, radii which result from the expansion in 
the SHO basis almost exactly reproduce values of finite element discretization in coordinate 
space. All calculations with the SHO basis have been performed with 20 oscillator shells. 
The same SHO basis has been used both for protons and for neutrons. 

In Fig. we display the proton and neutron density profiles for '^°Ne and ^°Ne. In our 
calculation ^°Ne is the last bound isotope. The proton density distributions practically 
do not change with the mass number, but the neutron densities display an abrupt change 
between ^°Ne and ^^Ne. A long tail reveals the formation of a multi-particle halo. Again 
we notice that calculations in the HO basis do not reproduce the asymptotic behavior of 
the neutron densities. Except for the radii shown for comparison in Fig. |T], we have used 
a HO basis of 20 shells. While the results of HO calculations display a typical Gaussian 
dependence at large distances from the center of the nucleus, SHO results with the same 
number of oscillator shells nicely reproduce the coordinate space calculations. More details 
are shown in Figs. |^ and ^, where neutron density profiles are plotted for the even isotopes 
30-40]\jg While all three methods (coordinate space, HO and SHO bases) produce almost 
identical density distributions for neutrons in the "inner" region r < 10 fm (Fig. results 
are very different for the asymptotic behavior (Fig. H). Not only is the SHO basis obviously 
superior when compared with results of the HO calculations, but to some extent it could 
have also an advantage over calculations in coordinate space. Since these are performed 
in a box of finite dimension (20 fm in the present calculations), near the boundary of the 
box a typical, very steep decrease in the densities is observed. Coordinate space solutions 
in this region of course have nothing to do with physical densities, rather they are an 
artifact of the finite dimension of the mesh. SHO solutions do not suffer from this drawback, 
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and the corresponding asymptotic densities are more natural. Of course one can increase 
the coordinate mesh, but this might not always be feasible, especially in more than one 
dimension. 

Calculations in the configurational space generally depend on the choice of parameters 
that determine the basis functions. In non relativistic self-consistent models the standard 
procedure is to minimize the energy of the ground-state with respect to the basis parame- 
ters [l^. The minimization is performed in each step of the self-consistent iteration. This 
method however cannot be applied in relativistic models. Due to the presence of the Dirac 
sea, the energy functional does not display a minimum, but rather a saddle point in the 
ground state. The usual approach then is to find regions of values of basis parameters for 
which all calculated quantities saturate, i.e. properties of the ground state become indepen- 
dent on the specific choice of the basis. In particular, the spherical HO basis is characterized 



by the oscillator length parameter b = yfi/muj. Calculations along the /5-stability line have 
shown that nuclear properties saturate for hcu ^ 41^4^/^ |[19[. The situation is different for 



nuclei at neutron drip lines. Properties associated with the halo states (as e.g. neutron 
densities and rms radii) depend linearly on the oscillator length parameter. The new SHO 
basis is determined by the two parameters: oscillator length b and the matching radius a. 
Our calculations have shown that quantities which characterize the ground-state saturate 
for those values of a and b for which the neutron densites decrease exponentially outside the 
nuclear volume. 

In Ref. [0 we have shown that the microscopic origin of the neutron halo can be found in 
a delicate balance between the self-consistent mean-field and pairing correlations. The RHB 
ground-state wave function can be written either in the quasiparticle basis as a product 
of independent quasi-particle states, or in the canonical basis as a highly correlated BCS- 
state. The canonical states are eigenstates of the RHB density matrix, and the eigenvalues 
are the corresponding occupation numbers. Since the density matrices in RHB are always 
localized, it follows that single-particle wave functions in the canonical basis are localized. 
The formation of the neutron halo is related to the quasi-degeneracy of the triplet of states 
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lf7/2, 2p3/2 and 2pi/2. The pairing interaction promotes neutrons from the lf7/2 orbital 
to the 2p levels. Due to their small centrifugal barrier, the 2p3/2 and 2pi/2 orbitals form 
the halo. In the present work we have also verified that calculations in the configurational 
representation with the SHO basis reproduce the energies and occupation probabilities of 
canonical neutron states calculated on the mesh in coordinate space. 

IV. SUMMARY 

In the present work we have applied the relativistic Hartree-Bogoliubov theory in the 
description of ground-state properties of neutron-rich Ne isotopes. The model has been ap- 
plied in the self- consistent mean-field approximation. Pairing correlations and the coupling 
to particle continuum states are described by the pairing part of the finite range Gogny in- 
teraction. In many applications of HFB and RHB it has been emphasized that, for a correct 
description of the asymptotic behavior of nucleonic densities, coordinate space solutions for 
the equations that describe the coupling between bound and continuum states are essen- 
tial. This is especially important on the neutron-rich side of the valley of /9-stability, where 
extremely diffuse and spatially extended neutron densities are found at the drip lines. For 
these nuclei solutions based on the expansion of wave functions and mean-field potentials 
in harmonic oscillator bases result in ground-state densities with wrong asymptotic proper- 
ties. Coordinate space calculations are relatively simple in one dimension, i.e. for spherical 
nuclei. In more than one dimension, numerical solutions on a coordinate mesh can become 
extremely involved, especially if finite range two-body interactions are considered. In fu- 
ture applications of the RHB theory to deformed exotic nuclei it would therefore be useful 
to have alternative methods for the solution of the coupled system of partial differential 
Dirac-Hartree-Bogohubov and Klein-Gordon equations. 

We have investigated properties of a new spherical basis in configuration space. This 
scaled harmonic oscillator basis is derived by a local scaling-point transformation of the 
spherical harmonic oscillator radial functions. The unitary scaling transformation produces a 
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basis with improved asymptotic behavior, while it preserves useful properties of the HO wave 
functions. The SHO basis has been employed in RHB calculations of ground-state properties 
of Ne isotopes. Results obtained in configurational representation have been compared 
with rms radii and densities calculated by finite element discretization in the coordinate 
space. It has been shown that densities calculated in the SHO basis nicely reproduce results 
obtained on the coordinate mesh. In particular, using a rather small number of oscillator 
shells, we have been able to reproduce the sharp increase in neutron rms radii which had 
been interpreted as evidence for the formation of a multi-particle halo. Neutron densities 
calculated in the SHO basis display the correct exponential decrease in the asymptotic 
region at large distances from the center of the potential. The present investigation has 
demonstrated that SHO bases can be used in the description of the structure of spatially 
extended nuclei. This might be especially important in studies of deformed halos, where 
coordinate space discretization algorithms become technically very complicated. 
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Figure Captions 

FIG. 1. Calculated proton and neutron rms radii for the Ne isotopes. Values obtained by finite 
element discretization on the coordinate mesh are compared with calculations in the configurational 
representation: expansion in the harmonic oscillator (HO) and scaled harmonic oscillator (SHO) 
bases. 

FIG. 2. Self-consistent RHB proton and neutron densities for the ground-states of ^°Ne and 
^°Ne. Density profiles calculated in configurational space are compared with results on the coordi- 
nate mesh. 

FIG. 3. Self-consistent neutron densities for the ^°~^°Ne isotopes. Distributions calculated by 
expansions in the HO and SHO bases are displayed together with profiles obtained on a mesh in 
coordinate space. 

FIG. 4. Same as in Fig. 3, but in logarithmic scale. 
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